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Abstract 

We assume that to zero order in the Knudsen number the deviation 
of the entropy from a background value vansishes. We then show that 
adding a super-Burnett term we obtain a stable state of rest. The resulting 
equations have the same form as the Burnett equations but with the value 
of some coefficients changed. In particular the result applies to nonlinear 
acoustics. 

We consider a slightly rarefied gas. To first order in the Knudsen number, 
Kn, the Navier-Stokes equations are valid. Burnett derived the corresponding 
equations to second order in Kn. Bobylev |5] showed that the state of rest is 
unstable for the Burnett equations, see also Uribe et al. [H]. In this contribution 
we make the assumption that the deviation of entropy from a background value 
is of the order of KnMa, where Ma, is the Mach number. This is the case for 
nonlinear acoustics, where Ma ~ Kn. We show that with an error Kn 3 the 
Burnett equations in this case can be replace by equations which are linearly 
stable. 

In the one-dimensional case, the Burnett expressions for the xx component 
of the pressure tensor P and the heat current q are, see Chapman & Cowling 
0] (the dots indicate nonlinear Burnett terms) 
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q = -kT x - - — (6 2 - 4 )v xx ... 
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We now linearize around a state at rest and uniform temperature and density, 
writing 

T = T [l + f], p = Po [l + p], v = J***±v. 

V m 

We introduce dimensionless variables, where the unit of length is of the order 
of the mean free path 
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X = X *«L I ™ t = t *^ m 
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In the sequel stars and tildes are omitted. We obtain the linearized one- 
dimensional Burnett equations 

Pt+v x = 0, (1) 
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V t = -(p + T) x + -V xx + -U! 2 pxxx - ^(^3 - U2)T XXX) (2) 
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-T t = -v x + -fT xx - -(0 4 - 9 2 )v xxx . (3) 

/ = 2mK/3fcfl/i is the Eucken number. In the calculations we use the value 
/ = 5/2. This is the lowest approximation in terms of Sonine polynomial 
expansion for any interatomic potential and is experimentally found to be a 
good approximation, see @j. 

Now we assume that the entropy change is of O(KnMa). We then have 

^ = ( 7 -l)^+0(XnMa). 
T p 

(7 = Cp/c v ). Linearizing and using dimensionless units we find 

(7 - l)p xxx - T xxx = O(KnMa). 

Hence, to within terms 0(Kn) we have for any value a which is 0(1) 
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u?.p xxx + (w 2 - uj 3 )T xxx = [ui2 + "(7 - l )]pxxx + g (^2 ~ w 3 - a)T xxx + O(KnMa). 

Thus, we can change the values of u> 2 , to 3 to u 2 och uj 3 in the linear part of the 
Burnett contribution. 

lj-2 = W2+a(7 — 1), (4) 
LU3 = 0J3 + cry. (5) 

Let us now choose a so that the coefficient of p xxx vanishes, or uj 2 — 0, This 
gives for a monatomic gas 
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LU 2 = 0, LJ 3 — UJ 3 - -W 2 -- 

For Maxwell molcules 

uj 2 = 2, uj 3 — 3; u 2 = 0, lj 3 = -2. 

For hard spheres 

lo 2 = 2.028, w 3 = 2.418, ub 2 = 0, uj 3 = -2.652. 
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Figure 1: Complex growth factor A for hard spheres, < k < 6. Rings Burnett, 
crosses our equations 



As a consequence, the p xxx term disappears. The sign of the T xxx term changes. 
For solutions proportional to exp[ikx + At] we find 
Asymptotically, for k — > oo we have, when (0 2 — 64) (u)2 — W3) > 0, 



A = 
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±n/^(^ - e 4 )(^ 2 - ^ 3 )fc 3 - (3/ „ +4) fc 2 
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Clearly, there is one mode that is nonpropagating and damped and there are two 
propagating, damped modes. One entropy mode and two sound wave modes. 
It is really not necessary to have 0J2 = 0, but just to have 0J2 — W3 > 0. 

Let us write down the resulting equations, first in the one-dimensional case. 
We neglect the nonlinear Burnett terms. 

Pt + {pv) x = 

p(v t + vv x ) = - — {pT) x + ^(pv x ) x + ^(u 3 - %u>2)^T xxx 
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^-p(T t + vT x ) = -- P Tv x + {kT x ) x + \{64 - 6 2 )—v xxx 
2m m dp 

Here, the coefficients of the Burnett terms can be taken at the background 
value, but the variations of p, k in the Navier-Stokes terms have to be taken 
inte account. 
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Figure 2: Phase veloctiy for < k < 1. Hard spheres. Diamonds Navier-Stokes, 
crosses our equations 

We now give the phase velocity. Note that the phase velocity is constant plus 
a term to order Kn 2 . Hence the deviation from straight lines for the Navier- 
Stokes equations is not physically relevant but that the deviation for the Burnett 
equations (and our equations) is. Note that the difference between the Burnett 
equations and our equations is for larger k than those shown in Fig. 2. 

The three-dimensional equations are 



In earlier contributions by Jin and Slemrod, [H] and by the present author 
jH], see also [B] and the Burnett equations were regularized to a set of 13 first 
order equations generally valid. 

The present regularization applies when the deviations of entropy is 0(KnMa) 
but gives equations which can more easily be applied for small Knudsen num- 
bers. The condition on the entropy applies for nonlinear sound propagation. 
The same assumption that Ma ~ Kn is called the weakly nonlinear case in Sonc 




p t + V • (pv) x = 0, 
p(v t + (vV)v) = --V(pT), 

TO 

+V • MVv + (Vv) T - ~(V ■ v)l]} + ^ 3 - ^ 2 )^A(VT) 
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v • VT) = pT(V • v) + V-(kVT) + -(0 4 - 2 )— A(V • v). 

m ip 
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UJ, where stationarity is assumed, but here sound waves are included as well. 

Recently, the present author has also obtained another set of regularized 
equations ^U], which like the Burnett equations and the equations in this work 
are equations for p, v, T . They are, however generally valid, with no limitation 
on entropy or Mach number. They are first order in time and third order in 
space, but also contain mixed derivatives first order in time and up to second 
order in space. 
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